Proteasome Inhibition Reprograms Chromatin Landscape in Breast Cancer.

The 26S proteasome is the major protein degradation machinery in cells. Cancer cells use the proteasome to modulate gene expression networks that promote tumor growth. Proteasome inhibitors have emerged as effective cancer therapeutics, but how they work mechanistically remains unclear. Here, using integrative genomic analysis, we discovered unexpected reprogramming of the chromatin landscape and RNA polymerase II (RNAPII) transcription initiation in breast cancer cells treated with the proteasome inhibitor MG132. The cells acquired dynamic changes in chromatin accessibility at specific genomic loci termed differentially open chromatin regions (DOCR). DOCRs with decreased accessibility were promoter proximal and exhibited unique chromatin architecture associated with divergent RNAPII transcription. Conversely, DOCRs with increased accessibility were primarily distal to transcription start sites and enriched in oncogenic superenhancers predominantly accessible in non-basal breast tumor subtypes. These findings describe the mechanisms by which the proteasome modulates the expression of gene networks intrinsic to breast cancer biology. SIGNIFICANCE: Our study provides a strong basis for understanding the mechanisms by which proteasome inhibitors exert anticancer effects. We find open chromatin regions that change during proteasome inhibition, are typically accessible in non-basal breast cancers.

APOBEC2 safeguards skeletal muscle cell fate through binding chromatin and regulating transcription of non-muscle genes during myoblast differentiation.

The apolipoprotein B messenger RNA editing enzyme, catalytic polypeptide (APOBEC) family is composed of nucleic acid editors with roles ranging from antibody diversification to RNA editing. APOBEC2, a member of this family with an evolutionarily conserved nucleic acid-binding cytidine deaminase domain, has neither an established substrate nor function. Using a cellular model of muscle differentiation where APOBEC2 is inducibly expressed, we confirmed that APOBEC2 does not have the attributed molecular functions of the APOBEC family, such as RNA editing, DNA demethylation, and DNA mutation. Instead, we found that during muscle differentiation APOBEC2 occupied a specific motif within promoter regions; its removal from those regions resulted in transcriptional changes. Mechanistically, these changes reflect the direct interaction of APOBEC2 with histone deacetylase (HDAC) transcriptional corepressor complexes. We also found that APOBEC2 could bind DNA directly, in a sequence-specific fashion, suggesting that it functions as a recruiter of HDAC to specific genes whose promoters it occupies. These genes are normally suppressed during muscle cell differentiation, and their suppression may contribute to the safeguarding of muscle cell fate. Altogether, our results reveal a unique role for APOBEC2 within the APOBEC family.

Reprogramming Chromosome Ends by Functional Histone Acetylation.

Cancers harness embryonic programs to evade aging and promote survival. Normally, sequences at chromosome ends called telomeres shorten with cell division, serving as a countdown clock to limit cell replication. Therefore, a crucial aspect of cancerous transformation is avoiding replicative aging by activation of telomere repair programs. Mouse embryonic stem cells (mESCs) activate a transient expression of the gene Zscan4, which correlates with chromatin de-condensation and telomere extension. Head and neck squamous cell carcinoma (HNSCC) cancers reactivate ZSCAN4, which in turn regulates the phenotype of cancer stem cells (CSCs). Our study reveals a new role for human ZSCAN4 in facilitating functional histone H3 acetylation at telomere chromatin. Next-generation sequencing indicates ZSCAN4 enrichment at telomere chromatin. These changes correlate with ZSCAN4-induced histone H3 acetylation and telomere elongation, while CRISPR/Cas9 knockout of ZSCAN4 leads to reduced H3 acetylation and telomere shortening. Our study elucidates the intricate involvement of ZSCAN4 and its significant contribution to telomere chromatin remodeling. These findings suggest that ZSCAN4 induction serves as a novel link between 'stemness' and telomere maintenance. Targeting ZSCAN4 may offer new therapeutic approaches to effectively limit or enhance the replicative lifespan of stem cells and cancer cells.

Chromosome-scale assembly of the streamlined picoeukaryote Picochlorum sp. SENEW3 genome reveals Rabl-like chromatin structure and potential for C4 photosynthesis.

Genome sequencing and assembly of the photosynthetic picoeukaryotic Picochlorum sp. SENEW3 revealed a compact genome with a reduced gene set, few repetitive sequences, and an organized Rabl-like chromatin structure. Hi-C chromosome conformation capture revealed evidence of possible chromosomal translocations, as well as putative centromere locations. Maintenance of a relatively few selenoproteins, as compared to similarly sized marine picoprasinophytes Mamiellales, and broad halotolerance compared to others in Trebouxiophyceae, suggests evolutionary adaptation to variable salinity environments. Such adaptation may have driven size and genome minimization and have been enabled by the retention of a high number of membrane transporters. Identification of required pathway genes for both CAM and C4 photosynthetic carbon fixation, known to exist in the marine mamiellale pico-prasinophytes and seaweed Ulva, but few other chlorophyte species, further highlights the unique adaptations of this robust alga. This high-quality assembly provides a significant advance in the resources available for genomic investigations of this and other photosynthetic picoeukaryotes.

Chromatin looping links gene expression to the assembly of transcription factories (Review).

Genes are not randomly dispersed within the nuclear space, instead they occupy precise sites either with respect to the nuclear lamina as well as to each other. This observation stands at the basis of the today well accepted concept of nuclear territories where any chromosome shows reproducible spatial connections with a selection of others in a general picture that meets a functional criterion where genes that answer the same stimuli are grouped in the same sites. In fact, transcription is not visible widely dispersed throughout the nucleus but is gathered in several 'granules', called transcription factories that accommodates ~10 genes concurrently transcribed. This dynamic behavior of chromosomes is allowed by changes in chromatin plasticity that are governed by several classes of proteins that either modify its building or induce post­translational modifications in the protein component of nucleosomes, triggering formation of chromosome loops that modify the location of specific sites along the DNA strand. For example, transcription associated to nuclear receptors benefits of the generation of nuclear ROS that induce nicks following activation of the DNA repair apparatus that enhance helix unfolding and chromosome bridging. In the present review, the role that protocols facing elucidation of chromosome architecture are playing and will play in the near future were highlighted in order to investigate composition of the transcription factories assembled in response of a specific trigger: The estrogen­sensitive transcription was cited but the authors are convinced that the same portrait will be observed with a multitude of (if not all) other stimuli.

Nuclear lamina component KAKU4 regulates chromatin states and transcriptional regulation in the Arabidopsis genome.

BACKGROUND: The nuclear lamina links the nuclear membrane to chromosomes and plays a crucial role in regulating chromatin states and gene expression. However, current knowledge of nuclear lamina in plants is limited compared to animals and humans. RESULTS: This study mainly focused on elucidating the mechanism through which the putative nuclear lamina component protein KAKU4 regulates chromatin states and gene expression in Arabidopsis leaves. Thus, we constructed a network using the association proteins of lamin-like proteins, revealing that KAKU4 is strongly associated with chromatin or epigenetic modifiers. Then, we conducted ChIP-seq technology to generate global epigenomic profiles of H3K4me3, H3K27me3, and H3K9me2 in Arabidopsis leaves for mutant (kaku4-2) and wild-type (WT) plants alongside RNA-seq method to generate gene expression profiles. The comprehensive chromatin state-based analyses indicate that the knockdown of KAKU4 has the strongest effect on H3K27me3, followed by H3K9me2, and the least impact on H3K4me3, leading to significant changes in chromatin states in the Arabidopsis genome. We discovered that the knockdown of the KAKU4 gene caused a transition between two types of repressive epigenetics marks, H3K9me2 and H3K27me3, in some specific PLAD regions. The combination analyses of epigenomic and transcriptomic data between the kaku4-2 mutant and WT suggested that KAKU4 may regulate key biological processes, such as programmed cell death and hormone signaling pathways, by affecting H3K27me3 modification in Arabidopsis leaves. CONCLUSIONS: In summary, our results indicated that KAKU4 is directly and/or indirectly associated with chromatin/epigenetic modifiers and demonstrated the essential roles of KAKU4 in regulating chromatin states, transcriptional regulation, and diverse biological processes in Arabidopsis.

Tissue-specific enhancer-gene maps from multimodal single-cell data identify causal disease alleles.

Translating genome-wide association study (GWAS) loci into causal variants and genes requires accurate cell-type-specific enhancer-gene maps from disease-relevant tissues. Building enhancer-gene maps is essential but challenging with current experimental methods in primary human tissues. Here we developed a nonparametric statistical method, SCENT (single-cell enhancer target gene mapping), that models association between enhancer chromatin accessibility and gene expression in single-cell or nucleus multimodal RNA sequencing and ATAC sequencing data. We applied SCENT to 9 multimodal datasets including >120,000 single cells or nuclei and created 23 cell-type-specific enhancer-gene maps. These maps were highly enriched for causal variants in expression quantitative loci and GWAS for 1,143 diseases and traits. We identified likely causal genes for both common and rare diseases and linked somatic mutation hotspots to target genes. We demonstrate that application of SCENT to multimodal data from disease-relevant human tissue enables the scalable construction of accurate cell-type-specific enhancer-gene maps, essential for defining noncoding variant function.

p53 rapidly restructures 3D chromatin organization to trigger a transcriptional response.

Activation of the p53 tumor suppressor triggers a transcriptional program to control cellular response to stress. However, the molecular mechanisms by which p53 controls gene transcription are not completely understood. Here, we uncover the critical role of spatio-temporal genome architecture in this process. We demonstrate that p53 drives direct and indirect changes in genome compartments, topologically associating domains, and DNA loops prior to one hour of its activation, which escort the p53 transcriptional program. Focusing on p53-bound enhancers, we report 340 genes directly regulated by p53 over a median distance of 116 kb, with 74% of these genes not previously identified. Finally, we showcase that p53 controls transcription of distal genes through newly formed and pre-existing enhancer-promoter loops in a cohesin dependent manner. Collectively, our findings demonstrate a previously unappreciated architectural role of p53 as regulator at distinct topological layers and provide a reliable set of new p53 direct target genes that may help designs of cancer therapies.

Study on chromatin regulation patterns of expression vectors in the PhiC31 integration site.

The PhiC31 integration system allows for targeted and efficient transgene integration and expression by recognizing pseudo attP sites in mammalian cells and integrating the exogenous genes into the open chromatin regions of active chromatin. In order to investigate the regulatory patterns of efficient gene expression in the open chromatin region of PhiC31 integration, this study utilized Ubiquitous Chromatin Opening Element (UCOE) and activating RNA (saRNA) to modulate the chromatin structure in the promoter region of the PhiC31 integration vector. The study analysed the effects of DNA methylation and nucleosome occupancy changes in the integrated promoter on gene expression levels. The results showed that for the OCT4 promoter with moderate CG density, DNA methylation had a smaller impact on expression compared to changes in nucleosome positioning near the transcription start site, which was crucial for enhancing downstream gene expression. On the other hand, for the SOX2 promoter with high CG density, increased methylation in the CpG island upstream of the transcription start site played a key role in affecting high expression, but the positioning and clustering of nucleosomes also had an important influence. In conclusion, analysing the DNA methylation patterns, nucleosome positioning, and quantity distribution of different promoters can determine whether the PhiC31 integration site possesses the potential to further enhance expression or overcome transgene silencing effects by utilizing chromatin regulatory elements.

Chromatin Immunoprecipitation to Investigate H2A.Z Dynamics in Response to Environmental Changes.

DNA methylation and posttranslational modifications of histones instruct gene expression in eukaryotes. Besides canonical histones, histone variants also play a critical role in transcriptional regulation. One of the best studied histone variants in plants is H2A.Z whose removal from gene bodies correlates with increased transcriptional activity. The eviction of H2A.Z is regulated by environmental cues such as increased ambient temperatures, and current models suggest that H2A.Z functions as a transcriptional buffer preventing environmentally responsive genes from undesired activation. To monitor temperature-dependent H2A.Z dynamics, chromatin immunoprecipitation (ChIP) of H2A.Z-occupied DNA can be performed. The following protocol describes a quick and easy ChIP approach to study in vivo H2A.Z occupancy.

The scales, mechanisms, and dynamics of the genome architecture.

Even when split into several chromosomes, DNA molecules that make up our genome are too long to fit into the cell nuclei unless massively folded. Such folding must accommodate the need for timely access to selected parts of the genome by transcription factors, RNA polymerases, and DNA replication machinery. Here, we review our current understanding of the genome folding inside the interphase nuclei. We consider the resulting genome architecture at three scales with a particular focus on the intermediate (meso) scale and summarize the insights gained from recent experimental observations and diverse computational models.

Integrative analysis of transcriptomic and epigenomic data reveals distinct patterns for developmental and housekeeping gene regulation.

BACKGROUND: Regulation of transcription is central to the emergence of new cell types during development, and it often involves activation of genes via proximal and distal regulatory regions. The activity of regulatory elements is determined by transcription factors (TFs) and epigenetic marks, but despite extensive mapping of such patterns, the extraction of regulatory principles remains challenging. RESULTS: Here we study differentially and similarly expressed genes along with their associated epigenomic profiles, chromatin accessibility and DNA methylation, during lineage specification at gastrulation in mice. Comparison of the three lineages allows us to identify genomic and epigenomic features that distinguish the two classes of genes. We show that differentially expressed genes are primarily regulated by distal elements, while similarly expressed genes are controlled by proximal housekeeping regulatory programs. Differentially expressed genes are relatively isolated within topologically associated domains, while similarly expressed genes tend to be located in gene clusters. Transcription of differentially expressed genes is associated with differentially open chromatin at distal elements including enhancers, while that of similarly expressed genes is associated with ubiquitously accessible chromatin at promoters. CONCLUSION: Based on these associations of (linearly) distal genes' transcription start sites (TSSs) and putative enhancers for developmental genes, our findings allow us to link putative enhancers to their target promoters and to infer lineage-specific repertoires of putative driver transcription factors, within which we define subgroups of pioneers and co-operators.

Early screening of colorectal cancer using feature engineering with artificial intelligence-enhanced analysis of nanoscale chromatin modifications.

Colonoscopy is accurate but inefficient for colorectal cancer (CRC) prevention due to the low (~ 7 to 8%) prevalence of target lesions, advanced adenomas. We leveraged rectal mucosa to identify patients who harbor CRC field carcinogenesis by evaluating chromatin 3D architecture. Supranucleosomal disordered chromatin chains (~ 5 to 20 nm, ~1 kbp) fold into chromatin packing domains (~ 100 to 200 nm, ~ 100 to 1000 kbp). In turn, the fractal-like conformation of DNA within chromatin domains and the folding of the genome into packing domains has been shown to influence multiple facets of gene transcription, including the transcriptional plasticity of cancer cells. We deployed an optical spectroscopic nanosensing technique, chromatin-sensitive partial wave spectroscopic microscopy (csPWS), to evaluate the packing density scaling D of the chromatin chain conformation within packing domains from rectal mucosa in 256 patients with varying degrees of progression to colorectal cancer. We found average packing scaling D of chromatin domains was elevated in tumor cells, histologically normal-appearing cells 4 cm proximal to the tumor, and histologically normal-appearing rectal mucosa compared to cells from control patients (p < 0.001). Nuclear D had a robust correlation with the model of 5-year risk of CRC with r2 = 0.94. Furthermore, rectal D was evaluated as a screening biomarker for patients with advanced adenomas presenting an AUC of 0.85 and 85% sensitivity and specificity. artificial intelligence-enhanced csPWS improved diagnostic performance with AUC = 0.90. Considering the low sensitivity of existing CRC tests, including liquid biopsies, to early-stage cancers our work highlights the potential of chromatin biomarkers of field carcinogenesis in detecting early, significant precancerous colon lesions.

MMGAT: a graph attention network framework for ATAC-seq motifs finding.

BACKGROUND: Motif finding in Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) data is essential to reveal the intricacies of transcription factor binding sites (TFBSs) and their pivotal roles in gene regulation. Deep learning technologies including convolutional neural networks (CNNs) and graph neural networks (GNNs), have achieved success in finding ATAC-seq motifs. However, CNN-based methods are limited by the fixed width of the convolutional kernel, which makes it difficult to find multiple transcription factor binding sites with different lengths. GNN-based methods has the limitation of using the edge weight information directly, makes it difficult to aggregate the neighboring nodes' information more efficiently when representing node embedding. RESULTS: To address this challenge, we developed a novel graph attention network framework named MMGAT, which employs an attention mechanism to adjust the attention coefficients among different nodes. And then MMGAT finds multiple ATAC-seq motifs based on the attention coefficients of sequence nodes and k-mer nodes as well as the coexisting probability of k-mers. Our approach achieved better performance on the human ATAC-seq datasets compared to existing tools, as evidenced the highest scores on the precision, recall, F1_score, ACC, AUC, and PRC metrics, as well as finding 389 higher quality motifs. To validate the performance of MMGAT in predicting TFBSs and finding motifs on more datasets, we enlarged the number of the human ATAC-seq datasets to 180 and newly integrated 80 mouse ATAC-seq datasets for multi-species experimental validation. Specifically on the mouse ATAC-seq dataset, MMGAT also achieved the highest scores on six metrics and found 356 higher-quality motifs. To facilitate researchers in utilizing MMGAT, we have also developed a user-friendly web server named MMGAT-S that hosts the MMGAT method and ATAC-seq motif finding results. CONCLUSIONS: The advanced methodology MMGAT provides a robust tool for finding ATAC-seq motifs, and the comprehensive server MMGAT-S makes a significant contribution to genomics research. The open-source code of MMGAT can be found at https://github.com/xiaotianr/MMGAT , and MMGAT-S is freely available at https://www.mmgraphws.com/MMGAT-S/ .

Two-way feedback between chromatin compaction and histone modification state explains Saccharomyces cerevisiae heterochromatin bistability.

Compact chromatin is closely linked with gene silencing in part by sterically masking access to promoters, inhibiting transcription factor binding and preventing polymerase from efficiently transcribing a gene. However, a broader hypothesis suggests that chromatin compaction can be both a cause and a consequence of the locus histone modification state, with a tight bidirectional interaction underpinning bistable transcriptional states. To rigorously test this hypothesis, we developed a mathematical model for the dynamics of the HMR locus in Saccharomyces cerevisiae, that incorporates activating histone modifications, silencing proteins, and a dynamic, acetylation-dependent, three-dimensional locus size. Chromatin compaction enhances silencer protein binding, which in turn feeds back to remove activating histone modifications, leading to further compaction. The bistable output of the model was in good agreement with prior quantitative data, including switching rates from expressed to silent states (and vice versa), and protein binding/histone modification levels within the locus. We then tested the model by predicting changes in switching rates as the genetic length of the locus was increased, which were then experimentally verified. Such bidirectional feedback between chromatin compaction and the histone modification state may be a widespread and important regulatory mechanism given the hallmarks of many heterochromatic regions: physical chromatin compaction and dimerizing (or multivalent) silencing proteins.

Transcription bodies regulate gene expression by sequestering CDK9.

The localization of transcriptional activity in specialized transcription bodies is a hallmark of gene expression in eukaryotic cells. It remains unclear, however, if and how transcription bodies affect gene expression. Here we disrupted the formation of two prominent endogenous transcription bodies that mark the onset of zygotic transcription in zebrafish embryos and analysed the effect on gene expression using enriched SLAM-seq and live-cell imaging. We find that the disruption of transcription bodies results in the misregulation of hundreds of genes. Here we focus on genes that are upregulated. These genes have accessible chromatin and are poised to be transcribed in the presence of the two transcription bodies, but they do not go into elongation. Live-cell imaging shows that disruption of the two large transcription bodies enables these poised genes to be transcribed in ectopic transcription bodies, suggesting that the large transcription bodies sequester a pause release factor. Supporting this hypothesis, we find that CDK9-the kinase that releases paused polymerase II-is highly enriched in the two large transcription bodies. Overexpression of CDK9 in wild-type embryos results in the formation of ectopic transcription bodies and thus phenocopies the removal of the two large transcription bodies. Taken together, our results show that transcription bodies regulate transcription by sequestering machinery, thereby preventing genes elsewhere in the nucleus from being transcribed.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans.

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

A fine-scale Arabidopsis chromatin landscape reveals chromatin conformation-associated transcriptional dynamics.

Plants, as sessile organisms, deploy transcriptional dynamics for adapting to extreme growth conditions such as cold stress. Emerging evidence suggests that chromatin architecture contributes to transcriptional regulation. However, the relationship between chromatin architectural dynamics and transcriptional reprogramming in response to cold stress remains unclear. Here, we apply a chemical-crosslinking assisted proximity capture (CAP-C) method to elucidate the fine-scale chromatin landscape, revealing chromatin interactions within gene bodies closely associated with RNA polymerase II (Pol II) densities across initiation, pausing, and termination sites. We observe dynamic changes in chromatin interactions alongside Pol II activity alterations during cold stress, suggesting local chromatin dynamics may regulate Pol II activity. Notably, cold stress does not affect large-scale chromatin conformations. We further identify a comprehensive promoter-promoter interaction (PPI) network across the genome, potentially facilitating co-regulation of gene expression in response to cold stress. Our study deepens the understanding of chromatin conformation-associated gene regulation in plant response to cold.

An integrated toolkit for human microglia functional genomics.

BACKGROUND: Microglia, the brain's resident immune cells, play vital roles in brain development, and disorders like Alzheimer's disease (AD). Human iPSC-derived microglia (iMG) provide a promising model to study these processes. However, existing iMG generation protocols face challenges, such as prolonged differentiation time, lack of detailed characterization, and limited gene function investigation via CRISPR-Cas9. METHODS: Our integrated toolkit for in-vitro microglia functional genomics optimizes iPSC differentiation into iMG through a streamlined two-step, 20-day process, producing iMG with a normal karyotype. We confirmed the iMG's authenticity and quality through single-cell RNA sequencing, chromatin accessibility profiles (ATAC-Seq), proteomics and functional tests. The toolkit also incorporates a drug-dependent CRISPR-ON/OFF system for temporally controlled gene expression. Further, we facilitate the use of multi-omic data by providing online searchable platform that compares new iMG profiles to human primary microglia: https://sherlab.shinyapps.io/IPSC-derived-Microglia/ . RESULTS: Our method generates iMG that closely align with human primary microglia in terms of transcriptomic, proteomic, and chromatin accessibility profiles. Functionally, these iMG exhibit Ca2 + transients, cytokine driven migration, immune responses to inflammatory signals, and active phagocytosis of CNS related substrates including synaptosomes, amyloid beta and myelin. Significantly, the toolkit facilitates repeated iMG harvesting, essential for large-scale experiments like CRISPR-Cas9 screens. The standalone ATAC-Seq profiles of our iMG closely resemble primary microglia, positioning them as ideal tools to study AD-associated single nucleotide variants (SNV) especially in the genome regulatory regions. CONCLUSIONS: Our advanced two-step protocol rapidly and efficiently produces authentic iMG. With features like the CRISPR-ON/OFF system and a comprehensive multi-omic data platform, our toolkit equips researchers for robust microglial functional genomic studies. By facilitating detailed SNV investigation and offering a sustainable cell harvest mechanism, the toolkit heralds significant progress in neurodegenerative disease drug research and therapeutic advancement.

The Myc-associated zinc finger protein epigenetically controls expression of interferon-γ-stimulated genes by recruiting STAT1 to chromatin.

The MYC-Associated Zinc Finger Protein (MAZ) plays important roles in chromatin organization and gene transcription regulation. Dysregulated expression of MAZ causes diseases, such as glioblastoma, breast cancer, prostate cancer, and liposarcoma. Previously, it has been reported that MAZ controls the proinflammatory response in colitis and colon cancer via STAT3 signaling, suggesting that MAZ is involved in regulating immunity-related pathways. However, the molecular mechanism underlying this regulation remains elusive. Here, we investigate the regulatory effect of MAZ on interferon-gamma (IFN-γ)-stimulated genes via STAT1, a protein that plays an essential role in immune responses to viral, fungal, and mycobacterial pathogens. We demonstrate that about 80% of occupied STAT1-binding sites colocalize with occupied MAZ-binding sites in HAP1/K562 cells after IFN-γ stimulation. MAZ depletion significantly reduces STAT1 binding in the genome. By analyzing genome-wide gene expression profiles in the RNA-Seq data, we show that MAZ depletion significantly suppresses a subset of the immune response genes, which include the IFN-stimulated genes IRF8 and Absent in Melanoma 2. Furthermore, we find that MAZ controls expression of the immunity-related genes by changing the epigenetic landscape in chromatin. Our study reveals an important role for MAZ in regulating immune-related gene expression.